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Knowledge  of  the  large-scale  vertical  velocity  is 
required  to  predict  the  evolution  of  the  atsospheric 
planetary  boundary  layer  (1PBL) .  since  naval  operations  are 
often  conducted  in  data  sparse  regions ,  single- station 
assesssents  of  the  vertical  velocity  are  necessary.  Data  to 
evaluate  such  assesssents  were  obtained  froa  ravinsondes 
taken  at  San  Kicolas  Island,  California.  Vertical  velocity 
estiaates  obtained  by  vertical  integration  of  the  aoisture 
budget  equation  (Q-aethod)  and  by  the  adiabatic  sethod  were 
used  to  initialize  an  1PBL  24-hour  prediction  aodel.  BBS 
error  statistics  on  predicted  inversion  height,  potential 
tenperature,  and  specific  huaidity  were  coaputed  for  the 
forecasts  and  coapared  to  the  BBS  errors  of  a  persistence 
forecast.  Calculation  of  the  vertical  velocity  by  the 


Q-aethod  shoved  the  aost  promise. 


However ,  no 


single-station  assessaent  sethod  improved  on  the  persistence 


forecasts. 
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I-  IHTBOD  DC  TIPS 


Cne  aspect  of  nasal  operations  is  battle  group  opera¬ 
tions  in  a  fixed  geographical  area.  This  is  currently 
referred  to  as  HODLOC.  Several  scenarios  lend  themselves  to 
HODLCC  operations:  the  Tiet-ffam  era  carrier  based  strike 

operations,  pre-positioning  for  a  projection  of  power,  or 
sisply  a  show  of  force.  Training  operations,  as  well,  are 
often  conducted  in  fixed  areas  to  avoid  hazards  to  coemer- 
cial  shipping  and  air  traffic  with  weapons  firings,  or  for 
privacy  in  operations. 

1  second  aspect  is  the  need  for  self  sufficiency  of  a 
deployed  force.  Any  antonoaons  ability  to  sake  environ- 
aental  prediction  within  the  battle  group  reduces  the 
requirements  for  shorebased  support,  decreases  delivery 
tiae  for  perishable  information  and  reduces  the  deaand  on 
coaaunication  assets. 

A  third  aspect  is  the  effect  of  the  atmosphere  on  the 
performance  of  weapons  and  sensors  that  depend  on  electro¬ 
magnetic  and  electro-optical  wave  propagation.  Heteorologi- 
cal  factors  which  affect  such  propagation  can  change 
significantly  in  four  to  six  hours.  The  importance  of  APBL 
prediction  to  naval  operations  lies  in: 
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•  Prediction  of  the  existence  and  evolution  of  electro- 
eagnetic  ducts. 

•  Prediction  of  the  inversion  height,  where  it  is  known 
that  optical  propagation  is  degraded  doe  to  turbulence. 

This  thesis  is  an  exaaination  of  the  accuracy  of  single¬ 
station  assessaents  of  the  large-scale  vertical  velocity, 
and  the  effectiveness  of  these  assessnents  for  an  existing 
&PBL  prediction  nodel.  Three  nethods  of  coaputing  the  ver¬ 
tical  velocity  will  be  discussed.  Vertical  velocities  cal¬ 
culated  by  single-station  nethods  will  be  used  in  an  &PBL 
aodel  to  ascertain  the  relative  aerits  of  these  nethods. 


II.  BACKS BOOH D 


The  single- station  forecasting  problea  probably  has  been 


of  interest  as  a  Military  requireaent  since  the  origin  of 


araed  conflict.  Oliver  and  Oliver  (1945)  pot  single-station 


assessaents  into  a  naval  perspective: 


Daring  the  last  fev  years,  vartiae  conditions  have  aade  it 
necessary  for  is  dated  coabat  units  to  issue  forecasts  in 
regions  where  no  network  of  a eteorological  stations  could 
be  available.  Frequently  the  data  froa  several  stations 
or  froa  reconnaissance  planes  are  available,  but  in  soae 
regions  the  forecaster  aust  rely  only  on  surface  and 
upper-air  observations  aade  at  his  own  station.  This  is 

Particularly  true  in  the  case  of  ships  at  sea.  Hence,  it 
s  important  to  develop  proficiency  at  extracting  inforaa- 
ticn  with  liaited  aerological  data....  He  aust  further  be 
able  to  distinguish  froa  solutions  which  are  internally 
consistent,  and  those  which  are  inconsistent  and  therefore 
to  be  discarded. 


In  iPBL  assessaents  the  assuaption  of  horizontal  hoao- 


geneity  (one-diaensionality)  is  often  used.  All  APBL  pre¬ 


diction  aodels  require,  as  an  externally  known  paraaeter. 


the  large-scale  vertical  velocity  which  direcly  affects 


changes  in  the  height  of  the  APBL.  Furtheraore,  the  height 


of  the  APBL  affects  other  properties  of  the  boundary  layer 


such  as  the  teaperature  and  hunidity. 


Data  for  this  exaaination  were  obtained  froa  ravinsondes 


taken  at  San  Hicolas  Island  (SKI),  California.  SKI  data  are 


sufficiently  representative  of  a  marine  environaent  for  this 
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purpose.  Battalino  at  al.  (1979)  have  atatad  that  condi¬ 
tion*  at  SHI  aay  ba  aarine,  continental,  or  aixad  depending 
on  tba  circulation  patterns,  k  decision  to  use  the  SHI  data 
vas  sada,  in  part,  because  of  the  availability  of  siailar 
data  froa  Vandenberg  (7 HD)  iir  force  Base,  California  and 
Point  Hugo  (PHTC) ,  California,  so  that  exaaination  of  the 
horizontal  variability  of  the  area  and  divergence  derived 
vertical  velocity  va.lues  could  be  accoaplished.  Pig.  1  is 
provided  for  geographic  orientation. 


Figure  1.  Southern  California_C  castal  Haters.  (After 
Battalino  et  al.  1979.) 
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A.  THE  PBEDICTIOH  MODEL 

In  APBL  prediction  eodel  in  which  the  aean  vertical 
notion  is  an  external  parameter  is  ased  to  establish  the 
affect  of  the  derived  values  of  vertical  velocity  on  boun¬ 
dary  layer  evolution.  The  aodel  used  is  based  on  the  one 
presented  by  Stage  and  Businger  (1981)  and  was  coded  for 
operational  use  by  Brower  (1982)  -  The  predictive  equations 
are: 

dq/dt  =  h**l(w Tqrr(0)  ♦  W.  Aq)  (la) 

d8./dt  *  h-‘  (w*0*  (0)  «•  Bc  -  Bb  ♦  B.A9.)  (1b) 

dh/dt  »  «D  ♦  ».  (1c) 

Where  q  is  the  total  water  content  (liquid  and  vapor)  ,  h  is 
the  inversion  height,  8.  is  the  equivalent  potential  temper¬ 
ature,  w *q ' (0)  is  the  surface  flux  of  moisture,  w*9£  (0)  is 
the  surface  buoyancy  flux,  A  is  the  difference  between  the 
variable  value  immediately  above  the  APBL  and  its  well  mixed 
value  in  the  APBL,  Bc  is  the  rate  of  radiative  heat  gain  per 
unit  area  at  the  cloud  base,  and  Rb  is  the  rate  of  radiative 
heat  loss  per  unit  area  near  the  cloud  top.  f,  is  the 
entrainment  velocity  and  W0  is  the  large-scale  vertical 
velocity  at  the  inversion  height. 
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The  eodel  ateoephere  is  two  layered,  coeposed  of  a  tur¬ 
bo  lent  lover  layer  in  vhich  equivalent  potential  tempera¬ 
ture,  specific  haaidity  and  vind  speed  are  constant  vith 
height,  and  a  stmiole  upper  ataosphere  vith  constant  specific 
huaidity  and  buoyancy  gradients.  The  layers  are  separated 
at  the  inversion  height  by  a  zero-order  discontinuity 
(juap).  Fig.  2  depicts  hov  a  saaple  sounding  is  character¬ 
ized  in  the  aodel  ataosphere. 

Surface  fluxes  of  aoisture  and  bouyancy  are  diagnosed 
froa  bulk  foraulae  (Lilly  1968).  The  liquid  vater  profile 
is  ccaputed  adiabatically  and  reduced  by  about  30  percent  to 
agree  vith  enpirical  results  (Fairall  et  al.  1981.)  Long 
vave  radiative  cooling  is  treated  using  the  Stephan-Boltzaan 
lav,  and  short  vave  radiative  varaing  is  treated  using  the 
delta-Eddington  approxiaat ion  (Fairall  et  al.  1981). 

The  large-scale  vertical  velocity,  VD  ,  is  treated  as  an 
external  parameter,  and  is  constant  over  the  eodel  run 
(2*-hr).  The  systea  is  closed  with  the  specification  of  the 
entrainaent  velocity  vhich  is  accoaplished  by  the  assuaption 
that  dissipation  of  turbulent  kinetic  energy  (TKZ)  is  a 
fixed  fraction  of  the  production  of  TKE  (Stage  and  Busin ger 
1981).  i  30- minute  tiae  step  is  used. 

in 
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Figure  2. 


Sasple  sounding  in  the  sodel  atsosphere.  Solid. 
Xi^e^n^resep-ts  data,  dashed  line  represents  the 


so 


itsosphere. 


III.  MET  HODS 


At  least  three  aethods  which  can  be  used  to  coapute  the 
large-scale  vertical  velocity  are: 

•  The  kineaatic  aethod. 

•  The  adiabatic  aethod. 

•  Integration  of  the  aoisture  budget  equation. 

The  kineaatic  aethod  requires  wind  sounding  data  froa  at 
least  three  appropriately  spaced  stations.  The  adiabatic 
aethcd  and  the  integration  of  the  aoisture  budget  equation 
are  truly  single-station  assessaents. 

A.  THE  KINEMATIC  METHOD 

e 

The  kineaatic  aethod  can  be  used  to  coapute  the  vertical 
velocity  when  nearly  siaultaneous  (within  two  hours)  sound¬ 
ings  are  available  at  three  or  aore  appropriately  spaced 
stations  such  as  VHD,  PHTC,  and  SNI  (Fig.  1).  Following 
Saucier  (1955),  the  vertical  velocity  can  be  derived: 

dw/dz  *  -  (du/dx  «•  dv/dy) ,  (2) 

where  u,  v  and  w  are  the  zonal,  aeridional  and  vertical  wind 
coapcnents.  The  horizontal  wind  coaponents  can  be  linearly 
apprcxiaated  by  a  Taylor  series  expansion  as: 


u  «  u.  ♦  (du/dx)  x  ♦  (du/dy)  y r 
▼  •  ▼.  ♦  (dv/dx)  x  ♦  (dv/dy)  y. 


(3a) 

(3b) 


Zg.  (3)  suggests  that  tha  horizontal  derivatives  aay  ba 
obtainad  froa  a  ra gras si on  of  tha  horizontal  wind  components 
onto  tha  horizontal  plana  (Fairall  at  al.  1981).  Decoaposi- 
tion  of  tha  wind  vector  froa  tha  sounding  data  at  tha  three 
stations  (SMI,  7 HD,  PHTC)  into  zonal  and  aeridional  coapo- 
nants  is  required.  Vertical  profiles  of  tha  horizontal  wind 
coapcnants  ware  constructed  at  a  200  eater  vertical  incre¬ 
ment  by  linear  interpolation  froa  the  irregularly  spaced 
data  levels  tc  those  levels  that  are  an  integer  multiple  of 
200  eaters.  Multiple  linear  regression  of  each  component  of 
the  horizontal  wind  onto  the  horizontal  plane  is  perforaed 
as  fellows: 


a  ♦  bx  ♦  cy. 

(«a) 

d  ♦  ex  ♦  fy. 

(4  b) 

where  a  through  f  are  the  regression  coefficients. 

Comparison  of  Eq.  (3)  with  Eq.  (S)  indicates  that  the 
coefficients  b  and  f  give  the  zonal  derivative  of  zonal  wind 
and  the  aeridional  derivative  of  aeridional  wind  at  the  cen¬ 
ter  of  a  triangle  foraed  by  the  three  locations  and  at  each 
vertical  level  where  the  regression  is  done. 

Vertical  integration  of  Eq.  (2)  is  perforaed  to  give  a 
vertical  profile  of  vertical  velocity  to  three  kilometres: 

3  km 

*<*)--  £<<*>)  ♦  (f)  )  Az  ,  (5) 
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where  w(z*0)  •  0  and  Az  •  200w 


B.  THE  ADIABATIC  METHOD 


For  isentropic  notions  there  is  no  tine  change  of  poten¬ 
tial  ten per a tare  following  the  notion.  Based  on  this 
assanption  Saucier  (1955)  specifies  the  adiabatic  set hod  for 
estinating  vertical  velocity  as: 


<ae/at)  ♦  v  (ve) 

(d8/a*> 


(6) 


where  9  is  the  potential  ten pera tare.  Saucier  further 
states  that  the  local  tenperature  change  due  to  advection 
nay  he  on  the  order  of  that  due  to  vertical  notion.  In  this 
study,  calculations  of  vertical  velocity  are  aade  both 
including  and  neglecting  the  adv active  tern.  Description  of 
the  advection  calculation  is  left  to  a  subseguent  section. 

The  vertical  gradient  is  deternined  by  a  linear  fit  to 
the  lapse  rate  of  potential  tenperature  above  the  inversion. 
The  local  tine  change  of  tenperature  is  obtained  at  fixed 
heights  (z  ■  1.0,  1.25,  1.5  kn)  froa  sequential  soundings  at 
the  sane  location.  These  heights  are  above  the  inversion. 
It  is  assuaed  that  the  vertical  velocity  at  the  surface  is 
zero.  The  vertical  velocity  at  the  inversion  is  taken  as 
the  average  of  the  vertical  velocity  at  the  fixed  heights 
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linearly  scaled  to  the  inversion  height.  Thus,  a  single¬ 
station  assessaent  of  the  vertical  velocity  is  obtained. 


C.  IRTEGBATIOH  OP  THE  HOISTOHE  BUDGET 

Integration  of  the  aoistnre  budget  agnation  to  cospnte 
the  vertical  velocity  was  performed  by  Lenschov  (1973). 
This  method  will  hereafter  be  referred  to  as  the  Q-aethod, 
and  is  based  on  the  assumption  of  veil-mixed  specific  humid¬ 
ity  in  the  boundary  layer,  and  that  changes  occur  due  to 
fluxes  at  the  sea  surface  and  inversion  only. 

Vertical  integration  of  the  aoisture  budget  equation 
through  the  boundary  layer  and  across  the  interface  between 
boundary  layer  and  the  stable  atmosphere  above  yields  the 
following  expression  for  vertical  velocity,  v,  at  the  level 
of  the  inversion,  h: 

4(9+1)  (dh/dt)  -  h(d(g*l)/dt)  ♦  w'q»  (0) 

w  »  ■■  . . .  -  . . —  .  1  — ,  (7) 

4(9+1)  ♦  C/Jh/2) 

where  g  and  1  are  vapor  and  liquid  water  contents,  w'g*  is 
the  aoisture  flux,  0  is  the  vertical  gradient  of  (g+1)  in 
the  aixed  layer,  and  4(9+1)  is  the  difference  between  total 
aoisture  in  the  nixed  layer  and  total  aoisture  immediately 
above  the  aixed  layer.  With  no  liquid  water,  z  «h,  and 
neglecting  aoisture  advection  Eg.  (7)  is  simplified  to: 
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v  * 


AqCdh/dt)  -  h (dq/dt)  ♦  w»q«(0) 

Aq 


(8) 


Lilly  (1968)  assumes  that  (3  is  zero.  Lenschov  (1973)  found 
/Jh/2  to  b«  saall  and  certainly  inch  less  than  the  aoisture 
jnsp. 

Iq.  (8)  provides  another  independent  assessment  of  the 
vertical  velocity  at  the  inversion.  The  tiae  derivatives 
(dh/dt  and  dq/dt)  are  determined  from  sequential  soundings. 
The  aoisture  jump  and  inversion  height  are  taken  as  averages 
for  the  two  soundings,  and  the  surface  aoisture  flux  is  det¬ 
ermined  by  the  bulk  method,  i.e. , 


a'q '  (0)  -  ChV10(q.  -  q10)  (9) 

The  exchange  coefficienx,  Ch,  is  computed  from  the  wind 
speed  dependent  relation  of  Large  and  Pond  (1980) .  The  sur¬ 
face  specific  humidity,  q,,  is  determined  via  Tetens*  for* 
aula  for  the  vapor  pressure  (Buck  1981). 


D.  DETERMINATION  OF  ADVSCTION  PBOH  THE  THERMAL  VIED 

Under  the  geostrophic  assumption,  the  horizontal  temper¬ 
ature  gradient  is  directly  proportional  to  the  vertical 
shear  of  geostrophic  wind.  Following  Byers  (1974): 

v  ■  f-»(d*/dx)p  :  a  *  #  (10) 


20 


where  v  and  u  ar«  the  meridional  and  zonal  components  of  the 
geostrophic  wind,  f  is  the  Coriolis  parameter,  *  is  the 
geopotential  and  the  horizontal  derivatives  are  on  the  con- 
stant  pressure  surface.  Vertical  differentiation  of 
Eq.  (10)  and  substitution  from  the  hydrostatic  equation 
yields: 

dr/dp  »  f~‘  (-a  )  ,  (11a) 

du/dp  ■  f-*M  a  )  ,  (11b) 

where  a  is  the  specific  volume.  Substitution  from  the 
equation  of  state  and  rearranging  terms,  the  components  of 
the  temperature  gradient  are: 

(dT/dx)p  -  -(fp/H  )  (dw/dp)  ,  (12a) 

(df/dy)p  -  (fp/H  )  (da/dp)  ,  (12b) 

where  B  is  the  gas  constant  for  dry  air. 

Bith  a  wind  sounding,  assumed  to  be  geostrophic  above 
the  boundary  layer,  the  horizontal  temperature  gradient  can 
be  specified  and  a  single- station  assessment  of  the  thermal 
advection  term  in  Eq.  (6)  is  obtained:  -  (V •  tfT)  *-  (V •  V9) ,  at 
a  constant  level. 

This  calculation  is  performed  by  subtracting  the  wind  at 
the  first  level  above  the  inversion  from  the  wind  at  an 
upper  level.  The  advect4“g  wind  is  the  mean  wind  in  the 
layer. 


E.  SIHGLB-STATIOM  1SSESSHBHT  SCOBES 

The  separate  vertical  velocities  coapated  froe  single- 
station  assessments  by  the  adiabatic  method  and  the  Q-method 
were  both  used  as  the  large-scale  vertical  velocity  in  the 
1PBL  model. 

Both  methods  required  use  of  finite  differences  in  time 
to  approximate  time  derivatives.  The  time  increment  used 
was  the  time  difference  between  subsequent  soundings  at  SHI; 
it  varied  from  four  to  eighteen  hours.  To  establish  the 
vertical  velocity  for  a  model  run,  the  time  weighted  average 
vertical  velocity  was  used,  as  computed  from  the  SHI  data 
for  the  24  hour  period  following  model  initialization. 

The  model  requires  initial  values  of  inversion  height, 
and  initial  profiles  of  potential  temperature  and  specific 
humidity.  Sea  surface  temperature,  required  for  computing 
surface  moisture  flux,  was  constant  through  the  model  day. 
Hind  speed,  required  for  computing  surface  moisture  and 
bouyancy  fluxes,  was  input  as  a  linearly  changing  function 
of  time. 

Three  model  runs  were  made  on  each  set  of  initial  condi¬ 
tions:  one  with  vertical  velocity  by  the  Q-method,  a  second 
with  vertical  velocity  by  the  adiabatic  method  and  the  last 
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with  vertical  velocity  set  to  zero.  The  advective  tera  of 
Eq.  (6)  it  was  found  to  be  aach  less  then  the  local  tiae 
change  of  potential  teaperature.  Thus,  it  vas  neglected. 

Boot  Bean  square  (BBS)  error  vas  coaputed  for  predicted 
inversion  height,  h,  potential  teaperature,  0,  and  specific 
huaidity,  q. 

BBS  error  «  1  Jd*'  “  >a 

where  B  is  the  nuaber  of  observations,  x,  is  the  observed 
value  of  a  variable,  and  y,  is  the  predicted  value  of  a 
variable.  Observed  values  were  obtained  froa  ravinsonde 

data  available  within  the  24  hours  following  the  tiae  of  the 
initial  data.  Typically,  two  or  three  soundings  were  avail¬ 
able  for  each  24- hour  period  following  initialization. 
Predicted  values  were  taken  froa  the  aodel  output  for  the 
tiae  co  rr.es  pan  ding  to  that  of  the  observed  values  (t  15 
ainutea,  due  to  the  30  ainute  aodel  tiae  step). 


1/2 


it.  U1QS11Z  sm&xiog 


Two  periods  (cases)  were  considered.  Case  I  was  froa  31 
lugust  to  4  September  1981  and  consisted  of  the  analyses  of 
nine  aodel  predictions  verified  against  twenty  observations. 
Case  II  was  froa  15  to  18  September  1981  and  consisted  of 
the  analyses  of  eight  predictions  verified  against  seventeen 
ob serrations. 

A.  CASE  I  (31  IDG  -  4  SEPT  1981) 

Throughout  the  period  the  eastern  North  Pacific  Ocean 
was  under  a  high  pressure  systea  that  was  moving  slowly 
eastward.  Southern  California  was  doainated  by  a  thermal 
trough*  leading  to  weak  offshore  flow  at  the  surface,  while 
California  was  to  the  south  of  the  strongest  westerlies  at 
upper  levels.  Pigs.  3  through  6  show  the  National  Meteoro¬ 
logical  Center  (NEC)  surface  and  500  aillibar  (ab)  analyses 
for  Case  I. 

Selected  prints  of  the  GOES  NEST  Infra-red  (IB)  iaagery 
indicate  low  broken  to  overcast  cloudiness  in  the  vicinity 
of  SVI  (Pig.  ?).  Surface  observations  at  SNI  indicated 
overcast  conditions  with  fog  at  night,  and  scattered  low 
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cloudiness  with  haze  persisting  after  the  late  aorning 
burn-off  of  the  fog. 

A  Pacific  high  pressure  systea  dominated  the  eastern 
ocean  area  west  of  California  on  31  August  (Pig.  3).  k  lov 
pressure  systea  was  located  in  the  Gulf  of  Alaska  with  its 
associated  weak  frontal  systea.  Tropical  Storn  Irwin  was 
dissipating  southwest  of  the  Baja  Peninsula.  A  theraal 
trough  doainated  the  central  California  valley.  very  weak 
westerly  flow  aloft  was  predoainant  over  southern  California 
and  aost  of  the  eastern  Pacific  Ocean. 

After  31  August  the  Pacific  high  intensified,  while  the 
low  pressure  systea  in  the  Gulf  of  Alaska  aoved  into  the 
northwest  United  States  and  the  assosiated  cold  front  dissi¬ 
pated.  The  theraal  trough  over  inland  California  aoved 
slightly  eastward  and  doainated  the  high  desert  region  of 
Arizcna  and  southeastern  California  (Pigs.  4  through  6). 
Flow  aloft  reaained  very  weak  over  southern  California  and 
the  low  level  flow  in  the  vicinity  of  SNI  was  very  weak  (3 
to  5  kt)  and  predoainant ly  offshore. 
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1249  GMT  31  AUG  81 


1315  GMT  1  SEP  81 


1246  GMT  2  SEP  81  1315  GMT  3  SEP  31 

Figure  7.  GOES  BEST  IB,  CASE  I. 
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B.  CASE  II  (15- 18  SEPT  1981) 

Throughout  the  period  the  polar  front  was  located  in  the 
eastern  part  of  the  North  Pacific  Ocean  and  Gulf  of  Alaska. 
A  high  pressure  systea  doeinated  the  area  southeast  of  the 
polar  front.  A  thermal  trough  was  located  over  the  inland 
California  valley  and  through  luch  of  the  southwestern 
United  States  and  Hexico.  Figs.  9  through  13  show  the  NSC 
surface  and  500  ah  analyses  for  Case  II. 

Satellite  imagery  indicates  considerable  clear  area  in 
the  vicinity  of  SHI.  Surface  observations  at  SHI  reported 
scattered  low  and  high  clouds  with  occasional  reduced  visi¬ 
bility  due  to  fog  and  haze.  Selected  prints  of  the  GOES 
WEST  IB  imagery  are  shown  in  Fig.  14  . 

Cn  14  September  high  pressure  prevailed  over  the  eastern 
Pacific  Ocean.  The  thermal  trough  dominated  the  southwest¬ 
ern  United  States  and  the  Hexican  Plateau.  A  low  pressure 
systea  with  a  weak  frontal  system  existed  in  the  Gulf  of 
Alaska.  Beak  unorganized  flow  existed  south  of  the  strong¬ 
est  westerlies  in  a  large  region  of  very  flat  height  gradi¬ 
ents  over  most  of  California  and  the  eastern  Pacific  Ocean 
(Fig-  9). 
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Ey  15  September  a  large  aaplitade  ridge  alo£t  bad  built 

( 

over  the  westers  Onited  States  and  Canada  (Fig.  10) .  The 
thersal  trough  weakened  in  central  California  as  the  ridge 
axis  in  the  eastern  Pacific  Ocean  Bigrated  eastward.  Low 
lewel  flow  in  southern  California  was  offshore  and  stronger 
(12  to  22  kt)  than  in  Case  I. 

C.  CXFFEBBVCSS  BETS SSI  CASES  I  AHO  II 

Conparison  of  the  aired  layer  and  overlying  data  clearly 
indicates  that  different  synoptic  situations  existed  for 
the  two  cases.  These  differences  are  sunaarized  in  Table  I 
and  illustrated  in  Fig.  8  : 


TABLE  1 

Coaparison  of  the  two  cases 


Case  I 

Case  II 

a  can  h 

604  a 

262  a 

range  of  h 

511  to  667  a 

145  to  357  a 

aean  6 

16.2  C 

18.1  C 

aean  A0 

15.0  C 

14.1  C 

aean  g 

8.3  g/kg 

9.8  g/kg 

aean  Aq 

-6.0  g/kg 

-4.4  g/kg 

aean  d8/dz  above 

h 

4.2  C/km 

2.8  C/ka 

aean  dq/dz  above 

h 

-.35  g/(kg-ka) 

-.83  g/(kg- 

On  the  average. 

inversion  height 

was  nearly  2. 

greater  in  case  I  than  in  case  II.  Case  II  had  a  waraer  and 
aore  aoist  boundary  layer  with  saaller  aagnitudes  of  juaps 


la  both  potential  taaparatara  and  specific  husidity.  Tha 
▼art leal  g radiant  of  potantial  taaparatura  (stability)  abova 
tha  invarsion  is  1.5  tiaas  graatar  ia  Casa  I  than  Casa  IX, 
sad  tha  vortical  soistora  gradient  is  naarly  2. 5  tiaas 
largar  in  Casa  I  than  Casa  II. 


Figure  8.  hr arage  aodel  Itsosphara:  Casa  I  vs.  Casa  II 
Solid:  Casa  I,  dash:  Casa  II 
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Figure  9.  Surface  and  500  mb  analyses,  0000  GMT  14  Sept  1981. 


Same  as  9  except  0000  GMT  18  Sept  1981. 


1015  GMT  15  SEP  81  0315  GMT  18  SEP  81 


Figure  14.  GOES  VEST  IR,  Case  II. 
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V-  8 ESQ ITS 


Direct  measurement  procedures  do  not  exist  to  obtain 
vertical  velocity  values.  Therefore  an  indirect  evaluation 
procedure  was  nesessary.  This  was  accomplished  by  using  the 
BBS  error  (section  III,E .)  for  aodel  predictions  of  inver¬ 
sion  height,  potential  tesperature  and  specific  humidity, 
with  the  large-scale  vertical  velocity  specified  by  the 
adiabatic  method,  Q-aethod,  vertical  velocity  equal  to  zero, 
and  a  persistence  forecast.  It  is  noted  that  the  Q-aethod 
is  based  on  aeasured  values  of  the  prediction  variables 
(Eq.  (1)).  Therefore  the  aodel  is  expected  to  perfora  well 
with  vertical  velocities  ccaputed  in  this  way. 

In  both  cases  the  inversion  persisted  for  several  days. 
Eq.  (1c)  illustrates  that  1  acts  to  keep  a  lid  on  the  boun¬ 
dary  layer  only  when  subsidence  (W  <0)  is  occurring  since  w 
provides  only  for  increasing  the  inversion  height,  i. e., 
"detrain sent"  is  not  allowed.  Then  froa  a  pragaatic  view, 
negative  vertical  velocities  (subsidence)  are  sought  while 
positive  ones  are  not.  Large-scale  subsidence  is  also 
expected,  frca  a  synoptic  viewpoint  in  both  cases,  since  SNI 
is  under  the  suaaertiae  subtropical  ridge. 


Q-method;  dash:  adiabatic  method;  double  dash:  kinematic 
method . 


Table  II  gives  the  averages  of  the  vertical  velocities 
for  each  case  calculated  by  the  Q-aethod,  the  adiabatic 
aethcd,  and  the  kineaatic  aethod.  The  standard  deviation  is 
given  in  parentheses.  Pigs.  15  and  16  shov  the  frequency  of 
occurrence  of  calculated  vertical  velocity  for  Case  I  and 
Case  II,  respectively.  Assuming  that  subsidence  existed  at 
the  inversion  height,  the  adiabatic  and  kineaatic  aethods 
are  unable  to  diagnose  subsidence.  The  Q-method  consis¬ 
tently  gives  subsidence  where  subsidence  is  expected.  Dis¬ 
cussion  of  possible  sources  of  errors  in  these  aethods  is 
left  to  a  subsequent  section. 

TABLE  II 

Vertical  velocity  (ca/sec)  comparisons 
Q-aethod  Adiabatic  method  Kineaatic  aethod 


Case 

I  -.26  (.14) 

0.0  (.22) 

1.2 

(1.7) 

Case 

n  -1.2  (1-1) 

.1  (-24) 

.7 

(1-3) 

A.  EISCUSSIOH  OF  CASE  I 

The  BBS  errors  for  Case  I  appear  in  Table  III  .  The 
average  tiae  from  initialization  to  verification  for  this 
case  was  15.6  hours.  Persistence  forecasts  yield  the  least 
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BUS  error  for  all  variables  (hf  9,  and  q)  .  When  the  aodel 
is  initialized  with  the  large-scale  vertical  velocity  com¬ 
peted  by  the  adiabatic  sethod  or  when  the  vertical  velocity 
is  asses ed  to  be  zero,  very  poor  results  are  obtained  in  the 
inversion  height  prediction.  Fertical  velocities  cosputed 
by  the  Q-aethod  result  in  better  prediction  of  inversion 
height  than  the  adiabatic  sethod. 

T1BLE  III 


BBS  Error  (Case  I) 


idia 


h 

Persistence  61s 

B«0  167  s 

Q-aethod  106  n 

oatic  sethod  249  a 


sethod  249  a 


9 

2.2  C 
2.2  C 
2.4  C 
2.2  C 


q 

.6  g/kg 
.8  g/kg 
.9  g/kg 
.7  g/kg 


The  large  errors  in  the  inversion  height  for  the  9*0 
prediction  occur  because  the  boundary  layer  is  constrained 
to  growth  by  entrainaent.  The  poor  results  of  the  adiabatic 
set bed  are  also  to  be  expected  because  ot  two  basic  probleas 
with  this  aethod.  First,  exaaination  of  Bq.  (6)  reveals 
that  the  lapse  rate  is  assumed  to  be  constant  and  the  numer¬ 
ator  is  approximated  by  a  finite  difference  in  tiae  over  six 
to  twelve  hours.  An  uncertainty  of  ±  one  degree  in  poten¬ 
tial  temperature  can  lead  to  an  uncertainty  in  vertical  vel¬ 
ocity  of  ±.5  (cm/sec),  assuming  nominal  values  of  the  free 

i 

ataoephere  lapse  rate  (2.5  to  5  C/Km) .  Since  the  actual 


44 


vertical  velocity  aay  be  positive  or  negative  and  often  on 
the  erder  .5  (ce/sec) ,  the  adiabatic  aethod  is  not  adequate 
to  specify  either  the  aagnitude  or  the  sign  of  the  vertical 
velocity.  Second,  one  aust  choose  an  appropriate  tine 
interval  over  which  the  finite  difference  applies.  In  prac¬ 
tice  the  tine  interval  vas  that  of  the  ravinsonde  schedule; 
however,  the  soundings  aust  be  scheduled  suffiently  close  in 
tiae  to  reveal  changes  in  sign  of  the  tendencies  of  the 
boundary  layer  variables. 

The  Q-aethod  also  suffers  froa  the  disadvantages  of 
finite  difference  in  the  specific  huaidity  (dq/dt)  and  the 
inversion  height  (dh/dt) •  However,  the  aoisture  flux  tera 
has  no  such  problea.  In  addition,  specification  of  the 
aoisture  and  inversion  height  variation  with  tiae  required 
for  these  calculations  does  not  require  profiling  devices 
such  as  radiosondes,  but  aay  be  obtained  froa  surface 
instruaentation. 

B.  EISCUSSICl!  Of  CASE  II 

lverage  tiae  froa  initialization  to  verification  was 
14.6  hours.  Table  IV  suaaarizes  the  results. 

lone  of  the  predictions  did  particularly  well  with 
respect  to  BUS  error  of  inversion  height.  The  sane  ordering 
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TABLE  I? 

BBS  Error  (Case  II) 


Persistence 

1*0 

.Q-aethod 
Adiabatic  method 


142  s 
527  ■ 


8 

1.5  C 
1.5  C 
1.0  C 
1.8  C 


.7gAg 
1 . 4g/kg 
2.4g/kg 
1.3#  A# 


of  BBS  errors  for  inversion  height  prediction  occurs  with 


the  vertical  velocity  specification  aethods  in  this  case. 


aajor  difference  in  these  results  is  that  the  model  pre¬ 
diction  of  specific  huaidity  with  vertical  velocity  by  the 
Q-aethod  exhibits  a  24ft  error,  such  greater  than  any  other. 
This  is  to  be  expected  since  the  average  inversion  height 
was  262a,  over  40ft  lover  than  in  Case  I,  and  in  Eg.  (la)  the 
specific  huaidity  tendency  is  inversely  proportional  to 
inversion  height.  Hence,  lov  initial  inversion  heights,  as 
vas  the  case,  coupled  with  predictions  of  inversion  heights, 
vhich  are  too  lov,  contribute  to  the  over  prediction  of  spe¬ 
cific  huaidity.  The  Q-aethod  yields  the  least  BBS  error  in 
potential  teaperature  predictions. 

C.  DISCUSSION  OP  THE  RIBEHATIC  BETHOD 

Vertical  velocities  calculated  by  the  kineaatic  net  hod 
vere  positive  aore  often  than  not  as  indicated  in  Pigs.  15 
and  16  .  They  vere  nuch  larger  in  aagnitude  than  those  cal¬ 
culated  by  either  the  adiabatic  aethod  or  the  Q-aethod. 


Therefore  application  of  these  to  the  aodel  would  have  pro¬ 
duced  larger  HAS  errors  than  any  presented.  Two  reasons  for 
the  poor  results  are: 

•  Iaproper  specification  of  the  bottoa  boundary  condi¬ 
tion. 

•  Snail  saaple  size  for  the  regression  approach. 

ever  the  ocean  the  bottoa  boundary  condition  for  Eg.  (2) 
is  that  the  vertical  velocity  equals  zero  at  the  sea  sur¬ 
face.  This  is  a  good  assuaption  since  the  surface  is  flat 
with  respect  to  a  local  earth  tanqent  plane.  However*  since 
SHI*  PHTC,  and  VHD  are  coastal  stations*  the  boundary  condi¬ 
tion  for  Eq.  (2)  requires  specification  of  the  terrain 
(Dutton  1976) . 

D.  COHCLOSIOH 

Since  both  Case  I  and  Case  II  represent  saall  sanples  of 
observations  the  error  statistics  have  little  stability. 
However*  they  do  suggest  that  vertical  velocity  calculated 
by  the  Q-aethod  has  the  aost  aerit  as  a  single-station 
assessaent.  clearly*  the  assuaption  of  no  vertical  notion 
or  vertical  notion  derived  fro  a  the  adiabatic  aethod  has 
questionable  value  in  1PBI  assessments.  Fig.  17  shows  a 
aodel  prediction  with  vertical  velocity  coaputed  by  the 
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Q-aethod.  Pig.  18  is  the  saae  model  run  but  with  upward 
vertical  velocity  coaputed  by  the  adiabatic  aethod.  Coapar- 
ison  of  these  figures  shows  the  good  agreeaent  in  q  pre¬ 
diction  indicated  in  Table  III  .  On  the  other  hand,  the 
superior  result  of  the  Q-aethod  is  well  illustrated  by  coa- 
parison  of  the  evolution  of  the  inversion  height.  The 
iapoxtance  of  diagnosing  subsidence  when  it  exists  is  evi¬ 
dent  by  this  coaparison. 

The  assuaption  that  SII  represents  a  marine  environaent 
is  a  possible  source  of  error  in  the  Q-aethod.  Also,  sea 
surface  teaperature  has  an  indirect  effect  on  the  calcula¬ 
tion  of  vertical  velocity  through  Eg.  (9).  The  sea  surface 
temperatures  used  were  obtained  from  the  Sea  Surface  Thermal 
Analyses  of  the  Hational  Oceanographic  and  Atmospheric 
Ada inist ration,  rather  than  being  in-situ  aeasureaents. 

Additional  research  is  needed  in  specification  of  the 
large-  scale  forcing  of  the  APBL .  Vertical  velocity  by  the 
Q-aethod  appears  proaising,  however  larger  samples  must  be 
obtained  to  gain  statistical  stability.  This  need  not  be  an 
expensive  venture.  One  only  needs  continuous  aeasureaents 
of  sea  surface  teaperature,  air  teaperature,  huaidity  and 
wind.  That  can  be  done  on  a  near-continuous  basis  with 


currently  available  instrumentation  that  does  not  place 
great  demands  on  personnel. 


o 
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Pigura  17. 
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